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Abstract 



(N 
> 

OO ' We present a measurement scheme for observing the Berry phase in a flux assisted Cooper pair 

pump — the Cooper pair sluice. In contrast to the recent experiments, in which the sluice was 
Tij- . employed to generate accurate current through a resistance, we consider a device in a supercon- 

^^ , ducting loop. This arrangement introduces a connection between the pumped current and the 

(^ I Berry phase accumulated during the adiabatic pumping cycles. From the adiabaticity criterion, 

we derive equations for the maximum pumped current and optimize the sluice accordingly. These 
results apply also to the high accuracy pumping which results in a potential candidate for a metro- 
logical current standard. For measuring the pumped current, an additional Josephson junction is 
^5 . installed into the superconducting loop. We show in detail that the switching of this system from 

Q I superconducting state into normal state as a consequence of an external current pulse through it 

may be employed to probe the pumped current. The experimental realization of our scheme would 
be the first observation of the Berry phase in superconducting circuits. 



* Electronic address: lmikko.mottonen@tkk.fil 



I. INTRODUCTION 

Adiabatic cyclic temporal evolution of any quantum system gives rise to geometric phases 
depending only on which path the system follows, not the speed of the cycle. For quantum 
states in degenerate subspaces, these geometric phases correspond to non-Abelian unitary 
operators within the subspace, also referred to as holonomies l|, |^. In holonomic quantum 
computation J3[ , these operators are of special interest since they represent quantum gates 
necessary for the actual computation. In the non-degenerate case, the geometric phase 
arising from a closed cycle in the adiabatic temporal evolution is referred to as Berry phase |^ 
corresponding to a phase of the quantum state under adiabatic temporal evolution. In 
general, the Berry phase is not directly observable and it is observed as a phase difference 
between two states which have traveled a different path during the adiabatic evolution, as was 
proposed in Refs. p, |6[ for the yet unrealized observation of the Berry phase superconducting 
quantum interference devices (SQUIDs). To date, the Berry phase has been observed |7|, for 
example, in systems of electrons [8| circulating about a wire carrying electric current known 
as Aharonov-Bohm effect, and in systems of neutrons p] or molecules [lO| circulating about 
a line of electric charge known as Aharonov-Casher effect. 

The peculiarity of the Berry phase in a phase biased array of Josephson junctions is that 
the geometric phase accumulated in the ground state of the Hamiltonian is closely related 
to the charge pumped thr oug h the device. This relation was first found by Aunola et al. 
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and also studied in Ref. 12|. Hence the yet unrealized observation of the Berry phase in 
superconducting circuits reduces in this structure into detection of weak currents 

We consider a flux assisted Cooper pair pump, also referred to as Cooper pair sluice [1 
consisting of a superconducting island separated by two SQUIDs which work as tunnel 
junctions with tunable Josephson energies. In order to achieve large and accurate pumped 
currents for a novel current standard, a voltage biased sluice was employed to pump current 
in Ref. |lj|. In contrast, we place the sluice into a superconducting loop, in which the device 
is phase biased, and hence may be employed to observe the Berry phase. 

For the measurement of the pumped current, we suggest that an additional Josephson 
junction is installed in the loop in parallel with the sluice. The pumped current modifies the 
switching statistics of the whole superconducting circuit into the normal conducting state 
under external current pulses. This change is shown to suffice for quantitative measurements 
of the pumped current. Furthermore, our measurement scenario resembles that of the su- 
perconducting qubit quantronium [l5|, |l6|, and hence we expect it to be experimentally 
feasible. 

Recently, a more sensitive measurement scheme than the one implemented in Saclay [l5|, 
16| for quantronium was implemented in Ref. [l7|. In this set up, the excitation and de- 



excitation of the additional tunnel junction introduces phase shifts to external current pulses, 
and hence no energy is transferred to the measured system. Whereas this method introduces 
improvements in the visibility and contrast of the measurement, its experimental implemen- 
tation is more demanding than that of the Saclay measurement scheme |l5Lll6l|. Moreover, 



we will show that the Saclay measurement scheme is actually more feasible in the observation 
of the Berry phase than in the quantronium measurement, and hence well justified. 

In the adiabatic temporal evolution of a quantum state, the parameters of the system 
Hamiltonian are controlled externally. The adiabaticity criterion states that the change in 
the parameters must be slow enough for the system to stay in the same instantaneous eigen- 
state of the Hamiltonian. This requirement limits the speed of pumping in the Cooper pair 



sluice, and hence the maximum pumped current. Thus the architecture of the sluice must 
be optimized in accordance with the limitations due to adiabaticity to obtain as strong a 
current signal as possible in the proposed measurement scheme. In Ref. jl3|, the errors in 
the pumped current due to the breaking of the adiabaticity were studied using a dynamic 
approach for a fixed design of the sluice. In contrast, we present both numerical and analyt- 
ical results for the maximum pumped current with any error rate. The results are obtained 
directly from the adiabaticity theorem and the spectral analysis of the system Hamiltonian. 
In addition, we apply the results to optimize the architecture of the sluice for both the Berry 
phase measurement and for high accuracy current pumping. 

The structure of this paper is the following: In Sec. m we present the Hamiltonian and 
the circuit of the Cooper pair sluice. Section IIIII is devoted to the study of the Berry phase 
and its connection to the pumped charge. Section HVl presents the optimization methods for 
the design of the sluice with respect to the adiabaticity criterion. The scheme to measure the 
pumped current is considered in Sec. El Section |VT1 summarizes and concludes the results 
of this paper. 

II. COOPER PAIR SLUICE 

The Cooper pair sluice consists of a superconducting island separated by two SQUIDs. 
Figure HI shows the sluice in a superconducting loop. In the case of vanishing loop inductance, 
the SQUIDs are equivalent to Josephson junctions with tunable Josephson energies £'ji($i) 
and £'j2($2); where $i and $2 denote the externally controllable magnetic fluxes through 
the first and second SQUID, respectively. On the other hand, the Coulomb energy for one 
excess Cooper pair to reside on the island is Eq = 2e^/Cs, where Cs represents the total 
capacitance of the island. If one neglects the parasitic capacitances, the total capacitance 
is given by Cs = 2Cj + Cg, where Cj is the capacitance of one of the identical SQUIDs and 
Cg is the gate capacitance which is used to charge the island with the gate voltage Vg. Thus 
the device can be considered to be a tunable Cooper pair transistor. 

The Hamiltonian of the sluice may be written in the eigenbasis of the phase operator 
as O 

H,i = Ec{n - UgY - Eni^i) cos(0 + y^/2) - E.^i^^) cos(0 - y^/2), (1) 

where h = —id^ is the number operator for the Cooper pairs on the island, rtg = VgCg/{2e) 
is the gate charge in units of 2e, (/? = y^i + ^92 is the phase difference over the device, and 
the phase on the island is given by = {ipi — ^2)12. The phase operators and y^j are 
related to the phase of the order parameter describing the many-body quantum state of the 
superconductor and should not be mixed with the Berry phase or the phase of the ground 
state of the above Hamiltonian. Since the sluice forms a superconducting loop, the phase if 
over the device is fixed by the magnetic flux $ through the loop as 

where $0 = vr/i/e is the flux quantum and we have assumed the the loop inductance is 
negligible. Thus ^p may be treated as a real number. 
The current operator of the fcth SQUID is defined as 

f 2ie^^ ^^^ 2edH 

4 = ^K,i/l = ^g^. (3) 




FIG. 1: Circuit for the Cooper pair sluice consisting of a superconducting island with gate ca- 
pacitance Cg, and two SQUIDs with magnetic fluxes <l>i and ^2 and phase differences of the order 
parameter tpi and ip2- The gate voltage is denoted by Vg. 



where fik = —id^^. is the Cooper pair number operator of the kth SQUID. For the sluice 
Hamiltonian Eq. (jSJ reduces into 



2eEjfc 

h 



sin((/9fc). 



(4) 



Equation (jH} ensures that for any states | ^(t)) and | x(^)) of the system the current operator 
possesses the property dt{^{t)\ (— 2e)rifc| x(t)) = (^(t)| j^l x(t)), i.e., the matrix elements of 
the current operator are given by the time derivative of the matrix elements of the charge 
operator. The critical current Ick is defined as the maximum current which can flow through 
the kth SQUID in the superconducting state. Hence, Eq. (JH) yields a relation 
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jfc 



'cfc 



$n 



(5) 



Thus no Cooper pairs can tunnel through the SQUID in the adiabatic evolution for van- 
ishing Josephson energy, and hence the SQUID is considered to be closed. In practice, the 
SQUIDs cannot be perfectly closed due to asymmetry in the tunnel junctions and finite 
loop inductance, which motivates us to define the minimum residual value of the Josephson 
energies of the SQUIDs to be E'^f^. In contrast, the maximum Josephson energy is denoted 
by Ey^^^. 

A typical pumping cycle of the sluice is described in Fig. |21 In the following, we assume 
that the SQUIDs can be perfectly closed, i.e., E^^ / E^^^ = 0. In the beginning of the cycle 
both of the SQUIDs are closed and the gate voltage is zero implying that the ground state 
of the Hamiltonian given by Eq. (Q) is an eigenstate of charge n with eigenvalue zero. Then 
the second SQUID is opened by increasing $2 to $o/2, while the first one is kept closed. The 
increase of the flux $2 is executed adiabatically such that the system stays in its ground 
state which is, however, not an eigenstate of charge h for finite $2- In the second step, 
the gate charge is brought to integer |25| n^^"^, after which the second SQUID is closed. 
At this middle point of the pumping cycle, the system is again in the charge eigenstate 
but corresponding to eigenvalue nf^^. Since the first SQUID has been closed for the whole 
evolution, we conclude that the n™*™ Cooper pairs have tunneled through the second SQUID 
to the island. Then we direct the n^^^ Cooper pairs through the first SQUID by opening 
it, taking rig to zero, and closing it. Thus n^^^ Cooper pairs have flowed through the sluice 
creating an effective current of Jp = 2e7if /T. Although the above arguments are valid in 



max g^j^j 



general, we will calculate analytically the pumped current in Sec. IIIII for Ec ^ E 

finite E^^. In analogy with a classical pump, the gate voltage in the sluice corresponds to 

a piston and the SQUIDs correspond to valves. 
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FIG. 2: Values of the externally controlled parameters £'ji($i), Ej2{^2), and rig of the Hamilto- 
nian given by Eq. (^ during an example pumping cycle of the Cooper pair sluice. The Josephson 
energies Eji and Ej2, and the gate charge Ug connected to the gate voltage V^ = 2eng/Cg are ex- 
pressed in terms of their maximum values E^^^ and n™'^^, respectively. The period of one pumping 
cycle is denoted by T. The different phases in the pumping cycle are separated by vertical lines. 
For e = EY^/EY^^ = 0, the pumping cycle is referred to as ideal, and hence either one of the 
SQUIDs is always perfectly closed. 



The maximum Josephson energy E^'^^ together with the adiabaticity criterion imposes 
upper bounds to the maximum pumped current, which we consider in detail in Sec. IIVI In 
the adiabatic temporal evolution, the SQUIDs stay in the superconducting state, and hence 
the average current through the sluice is given by 



2e dH,i 

h dip 



dipi d(p2 



h 



(6) 



The above equation yields a rather large upper bound for the effective pumped current, i.e.. 



the maximum critical current J™^^ 



27r^max/$^ of the SQUIDs. 



III. MANIFESTATION OF THE BERRY PHASE IN COOPER PAIR PUMPS 



Let us concentrate on the relation between the accumulated Berry phase and the pumped 
current. In the case of one- dimensional arrays of Josephson junctions, the expression of the 
pumped charge was first obtained in Ref. [18[ and its relation to the Berry phase was found 
in Ref. |ll|. Work concerning not only the Berry phase but also non-Abelian phases and 
their connection to the pumped charge in rather general Josephson devices is currently in 
preparation [19!]. 

In Sec. IIIIAl we consider a more general system than the sluice presented in Fig. Q Here 
we assume in excess of cyclic and adiabatic temporal evolution only that the average current 
operator / = 2e{d^pH) / h is well defined, i.e, that the Cooper pair pump has only one input 



and one output. We derive the connection between the pumped charge and the Berry phase 
for this rather general pump. In Sec. IIIIBl we focus on the sluice and calculate the pumped 
current using a two-state approximation which is valid in the limit Eq ^ Ef^^- 

A. Connection between the pumped charge and the Berry phase 

In adiabatic temporal evolution, an eigenstate of the system Hamiltonian H is slowly 
varied with respect to a set of real parameters {g/c}, here denoted by a vector q. The initial 
state of the system at t = to must be an eigenstate | m; q) of the Hamiltonian, for which one 
obtains 

^(q) I m; q) = e^(q) | m; q) , (7) 

where em(ct) is the eigenenergy of the corresponding state. In this section, we assume that 
the system is very accurately in its non-degenerate eigenstate |r;q). For these so-called 
zeroth order assumptions, the state of the system at any time instant is described by J20| 

|V^o(t)) = e^^^W|r;q), (8) 

where the phase 9r{t) = 6dr{t) + ^gr(^) is a sum of the dynamic phase 

Oar = -^ f\r; q(r) | H\ r; q(r)) dr = -^ l" e,(q(t)) dr, (9) 

and the geometric phase 

r-q(t) 



Ogr = i (r;q(r)|9^|r;q(r))dr = / (r; q| Vq| r; q) ■ dq. 
Jto Jciito) 



(10) 



In cyclic temporal evolution, the dynamic phase is proportional to the period T of the cycle, 
but the Berry phase, defined by 

6'Br = « ^(r;q| Vq|r;q) -dq, (11) 

only depends on the contour q = {q(t)|to < t < T -|- to}, according to which the state | r; q) 
is varied in the cycle. 

To obtain the connection between the Berry phase and the pumped current (l8l.ll9J. we 
consider the leading order corrections to the state vector | ^(t)) in the speed of the pumping 
cycle |q|. Hence, we express it as 

CO 

fc=0 

where Ck{t) are complex numbers. By differentiating Eq. (J12j) with respect to time and 
employing the Schrodinger equation we obtain 

cut) = - E Cfc(t)e*['''=W-^-Wl (m; q| d,\ k; q). (13) 
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Since we control our system such that it remains almost exactly in the state | r; q), we may 
neglect all other terms in Eq. (fT^ except the one with k = r. Thus we obtain for m ^ r 
the first order correction to the coefficients Ck as 

Cm{t) = - f {m; q| dt\ r; q)e^['''-(*')-^-(*')] dt'. (14) 

Jto 

As shown for example in Sec. IIVI the absolute value of the above coefficients Ck is propor- 
tional to |q| in adiabatic and cyclic temporal evolution. 

The total mean charge through the system in one pumping cycle may be written with 
the help of the average current operator / as 

Qtot= /" m)\i\i^{t))dt. (15) 

The contribution to the total charge purely from the zeroth order term presented in Eq. (jHJ 
is observed to originate from the usual supercurrent and is denoted by 



rto+T 
Qsr= (r;q(t)|J|r;q(t))dt. (16) 

J to 

ped charge Qp := Qtot — Qs assumes the form 
Qp, = 23fJe^*°'''' (5^C,(t)(r;q(t)|/| fc; q(t))e*[^'=W-^^WM dt, (17) 



-to+T 
'to 

On the other hand, the pumped charge Qp := Qtot — Qs assumes the form 

-to+T 
'*0 \kj^r 



where we have neglected the terms containing \Ck\'^ which vanish in the adiabatic limit. By 
inserting Eq. (fT^ into the above equation and changing the order of the resulting integrals 
we obtain 



fto+T 

,'\\„i[Or{t')-ek{t')] 



{{k-q{t')\d,\r-q{t'))e' 

j-to+T 

X / (r; q(t) | /| fc; q(t))e'[^^W-^'- «] dtjdt'. (18) 



In the above integral with respect to t, the integrand may be expressed as 

^ |zn(r; q(t)| /| k- q(t))e*[^^«-^^Wl/(^fc " ^.)} + 0{H), (19) 

The term 0(|q|) in the above equation may be dropped since it becomes negligible in the 
adiabatic limit compared with the first term. Hence the pumped charge is obtained to be 

Qp.= in^u^y \ r^^ ^""^ "^^^^l ^' ^' "^^^^^^'^ "^^^^l ^*' ""' "^^^^^ dt 

^ Ckjto + T){r; q(to + T)\i\ k; q(to + T))e^[^^fa+^)-^-fa+^)] i 

Er — Sk -I 



Since the excitation amplitudes Ck tend to zero in the adiabatic hmit, we obtain the known 
equation [18| for the pumped charge 






(21) 



which only depends on the contour g, according to which the state | r; q) traverses, not the 
speed |q| of the pumping cycle. Although the above equation contains an integral over a 
closed contour, it is also valid for open contours. Furthermore, the matrix elements of the 
average current operator / = 2e[9^, H]/h can be expressed as 



(r; q| /| k] q) = ''- -{r; q| d^\ k; q). 



h 



which yields together with Eq. (PT|) 



Qpr = -4eS5m (h (r; q| 9^Vq| r; q) ■ dq. 



(22) 



(23) 



In the above derivation, we have employed the fact that (r; q| d^\ r; q)(r; q| Vq| r; q) is real 
and that the basis {| fc; q)} is normahzed spaning the whole configuration space. 

Equation ()23p for the pumped charge resembles Eq. ()11|) except for the operator d^. This 
observation motivates us to study the derivative of the Berry phase with respect to the phase 
difference across the sluice (p which is fixed by the total magnetic flux $ as in Eq. Q, and 
hence may be treated as a scalar variable. Hence, Eq. pi|) yields 



dj 



ifi^Br 



dq. 



{d^\ r; q))^ Vq| r; q) + (r; q| (9<^Vq| r; q) 

The first term in the right side of Eq. ()24j) may be written in a form 

(a^|r;q))"^Vq|r;q) = Vq (5<^| r; q))"^ | r; q) - (Vq5^| r; q))^ | r; q). 



(24) 



(25) 



The integral of the term Vc 



{d^\r;q}y\r;q) 



theorem. Thus Eq. (j24p may be recast into 



over a closed contour q vanishes by the Stokes 



d,„9 



V^^Br 



-2Q=m ® (r; q| 9^Vq| r; q) ■ dq. 



(26) 



and hence the comparison of Eqs. (PHj) and (j^Uj) yields the relation between the Berry phase 
and the pumped charge 

Qp = 2e9^^B. (27) 

The interpretation of Qp is the average pumped charge through the device, since it arises 
from the average current operator I. Furthermore, if the pump consists of a linear chain 
of islands separated by SQUIDs, the charge conservation and cyclic temporal evolution 
assures that the pumped charge through all the SQUIDs separately equal to Qp |l8|. Since 
the Berry phase is not directly observable in a one- dimensional subspace of any quantum 
system, only the derivative of the Berry phase with respect to the total phase difference 
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appears in Eq. ()27p. We do not regard this property to be a limitation in measuring the 
Berry phase in Josephson circuits, but rather equate it with the fact that generally the Berry 
phase is detected as a phase difference between two orthogonal quantum states. In contrast 
to the Berry phase, the charge corresponding to the dynamic phase is related to the usual 
supercurrent through the device as 

Qs = -2ed^9d, (28) 

which is obtained by differentiating Eq. © with respect to if and comparing the result with 
Eq. (|T^. 

B. Berry phase in the Cooper pair sluice 

Assuming Ec ^ Ef'^'^ and Ug G [0, 1], the state of the Cooper pair sluice can be described 
by two eigenstates | 0) and 1 1) of the charge operator n corresponding to eigenvalues and 1, 
respectively. We employ this approximation to demonstrate the results of Sec. 1111 Al for the 
sluice up to linear terms in 

pmax 

5 := -j^. (29) 

In the discussion of Sec. Ull we found that the pumped charge is 2en^^^ in an ideal pumping 
cycle shown in Fig. |21for ^^es^^max _ g^ j^^ contrast, we consider here SQUIDs with finite 

— ^' (30) 



TTTnax 



We will, however, assume that the residual Josephson energy is small, and hence study its 
effect only up to linear terms in e. 

From the expression h = —ids of the number operator, we observe that 



e^"^\n) = \n±l), (31) 

where | n) is an eigenstate of the number operator h with an eigenvalue n. Thus we de- 
compose the trigonometric functions of Eq. (^ into exponential functions and obtain the 
approximate two-state Hamiltonian in the basis {| 0), | 1)} as 

H -( Ecnl _^e-./2 _ ^e^./2X 

^2s - ^_^gi^/2 _ ^e-^V^/2 E^^^_ ^J2 J ■ l-^^J 

The ground state | g) and the excited state | e) are obtained by diagonalizing this Hamilto- 
nian, and can be expressed as 

\g) = e'^a\0) + b\l), 

|e) = e^''6|0)-a|l), (33) 

where the amplitudes a and b are positive real numbers satisfying a"^ + b"^ = 1 and 7 is the 
relative phase difference of the charge states. The overall constant phase of the eigenstates 
is physically unobservable and can thus be chosen as in Eq. (jHHj) . The amplitudes a and b 
are obtained to be 

a^ = l-b^ = }L[l- ^ 1, (34) 



where we have defined rj = rig — 1/2 and £'12 = ^aZ-Eji + -E'j2 + 2EjiEj2 coscp. Moreover, 
the phase difference 7 assumes the form 

7 = arctan I — -— tan — ) . (35) 

Our aim is to compute the pumped charge utihzing Eq. ()21|) . and compare it with the 
derivative of the Berry phase with respect to the phase difference if. In this two-state 
approximation, the sum in Eq. pijl consists of only one term 



Qp2. = 2h^m i iSLSJj^^^e; q| VJ g; q) ■ dq, (36) 



? ^g Ee 



where we have employed the fact that all the current operators /, Ji, and I2 result in 
equal pumped charges due to the charge conservation. The symbols Eg and ^e denote the 
eigenvalues of the Hamiltonian H2S corresponding to the states | g; q) and | e; q), respectively. 
The parameters q are chosen to be q = (0^,7). With the help of Eqs. (jlj and (|^. the 
matrix element of the current operator I2 in Eq. (jHH|) is expressed as 

(g; q| /2I e; q) = jE,2[e'^''-^/%' + e-^^''-^/'^ a% (37) 

We write the gradient part in Eq. (jHUj) as 

(e;q|Vq|g;q)-dq = ^a6d7 + ^^, (38) 

and the energy difference as 

E,-Eg = 2Ec ^J Elr^-^ + El^. (39) 

By a substitution of Eqs. (jHTj) . (jSHl), and ^^ into Eq. (jHHjl we obtain 

Q,. - -2e*e^f|e.<--/^..^ + ,-(,-./.,.] I^Jj^f^ + id(a^)]. (40) 

Since the part containing d7 in the above integral vanishes as 0{5'^), it is neglected in this 
approximation. With the same argument we neglect the terms containing d(a^) for the 
parts of the contour <; in which E'ji or i5j2 is varied. Furthermore, we employ the identity 
6^ = 1 - a^ to Eq. (gOI) and obtain 

Qp2s = -2e I ^ cos(7 - v/2) d(a2) + (9(5^), (41) 

where the integration is over the second and fifth pumping phase, see Fig. El For the second 
pumping phase, ii^j2 = Ef^^^, and up to linear order in e we have E12 = EJ^''^^{1 + ecosip)/2 
and cos(7 — f/2) = 1. Furthermore, a^ varies from 1 — (9(5^) to 0(5^), and hence the 
contribution to the pumped charge from this domain is 2e(l — e cos if) . For the fifth pumping 
phase, Ej2 = Ef", E12 ^ Ef'^^{l + ecosip)/2, and cos{'j — ip/2) ^ cos(99). In this part of the 
cycle a^ varies from 0{6^) to 1 — 0{6'^), and hence the contribution to the pumped charge 
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up to linear order in e is given by — 2eecos(y9. Thus the pumped charge of the sluice in the 
adiabatic cycle of Fig. |21 assumes the form 

Qp2s = 2e (1 - 2e cos(/^) + 0{6^) + 0{e^). (42) 

For vanishing residual Josephson energy, the above equation reduces to the general result 
for ideal pumping cycle Qp = 2eri^^^, where the number of Cooper pairs pumped in one 
cycle n™*^^ is unity. 

To show the relation between the pumped charge and the Berry phase, we employ Eq. ()11|) 
and write 



du9.K — i 



7B2s = '^ i (g; q| Vql g; q) ■ dq = - i a^ d7, (43) 

where the vanishing integral involving the differential d(a^) is not shown. Since the param- 
eter 7 does not depend on the gate charge, the Berry phase is accumulated only when the 
Josephson energies are changed in the cycle of Fig. |21 Since a? = (9(5^) for rj = 1/2, the 
only non- negligible contribution to the integral in Eq. (j43p is obtained in the first and last 
pumping phase, for which a^ = 1 + 0(6"^). Thus the Berry phase can be expressed as 

eB2s = ^-2e sin((^) + 0{6^) + 0{e^). (44) 

Indeed, Eqs. fl^^ and (jii|) demonstrate that 

Qp2s = 2e9<^6'B2s, (45) 

which demonstrates for the sluice the validity of the general result given by Eq. (j27p . 

IV. ADIABATICITY CRITERION FOR THE COOPER PAIR SLUICE 

In the above sections, we have assumed that the sluice is more or less perfectly in its 
instantaneous ground state during the temporal evolution. However, this assumption is 
strictly speaking valid only if the changes in the parameters q of the system Hamiltonian 
are infinitely slow. On the other hand, we wish to repeat the pumping cycle as fast as possible 
for a given n^^^ to create a detectable current. Thus a quantitative study of the pumping 
errors due to finite pumping rate is required. In this section, we derive upper bounds for the 
population of the excited states of the sluice in the spirit of adiabatic theorem. Moreover, 
we apply the results to the case 5 <^ 1 for which the two-state approximation is vahd, and 
to the case S ^ 1 for which we employ a perturbed harmonic oscillator approach. In these 
two cases, we also derive the maximum pumped current for small pumping error rate. 

As in Sec. IIII Al we study a rather general quantum system described by the Hamilto- 
nian -ff(q), parameters q, and the instantaneous eigenbasis {\m;q)} defined in Eq. (jTj). We 
consider a representation of the state of the system given by Eq. (jl2|) . Since the aim of this 
section is to show how we can retain the system in its ground state, i.e., |Co(t)P ~ 1, we may 
employ Eq. (fT^ with r = to describe the time evolution of the excitation coefficients Ck- 
It follows from Eq. ((7j) that 

(m;q|gf|fc;q)= . (46) 
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By insertion of Eq. ()46|) into Eq. ()14|) we obtain 



'{m;q\ideH)\0;ci) 



cut) = / ^"^'^'7 ; ^^'''"^"^''"^"^' di'- (47) 

Let us assume that we adjust the speed of the temporal evolution such that 
(m; q| {dtH)\ 0;q} /[emit) — eoit)]"^ is constant in time. Furthermore, if we neglect the ge- 
ometric part of the phase and assume that the energy separation is time independent, we 
obtain an estimate 

'-'*" <- -^i3iar • 

Although the above assumptions may be considered to be crude in general, they turn out to 
be accurate for a sluice Hamiltonian with at least one of the Josephson energies much larger 
than the Coulomb energy Eq. 

Without any other assumption about the temporal evolution of the system except that 
it is cyclic and nearly adiabatic, Eq. (^Tj) yields 

1 - z/^+i 
CmiNT + to)=CmiT + to) "^ , forz/„^l, (49) 

where A^ is a positive integer and u^ = e*'^'"*^^-'"^'^^^-". Since the dynamic phase depends 
explicitly on the period T and the geometric phase only on the contour q, according to 
which the parameters q vary under the cyclic temporal evolution, the period T can be 
chosen for example such that Um = i for a fixed m, which implies \Cm{NT)\ < y2\Cm{T)\. 
Thus the probability for the sluice to be in the excited state | q; m) is bounded. The above 
discussion is, however, strictly valid only for the mth state, whereas u^ can generally be 
very close to unity for k ^ m. In the extreme case z/^ = 1, the amplitude of the kth state 
increases linearly in time, and hence the adiabaticity is lost if there is no dissipation in the 
system. To take the dissipation into account, we estimate from Eq. ()47p that 

^ < IC"^'»l(f.^)|0;'i)l < _ l<-^'^lf^)l°;-»>l :. .„. (50) 

emit) - eoit) 0<t<T Emit) - eo(^) 

For the excitation probability Pmit) = |Cm(^)P, Eq- ff^ and addition of relaxation to the 
system yields 

Pm < AKm^/Pmit) - r„P„(t), (51) 

where F^ is the relaxation rate of the state | q; m). Equation ^^ yields a bound iAKm/Tm)'^ 
for the excitation probability Pmit), and hence 

\CUt)\<P^. (52) 

From Eq. ()52|1 we can conclude that if the relaxation is dominant, i.e., F^ > 2(5^ — ^o)/^, 
bound in Eq. pH|l is valid even though the assumption that the energy difference Em — ^o is 
independent of time fails. Moreover, since k^, tends to decrease rapidly with increasing m, 
it may be argued that the excitation amplitudes of the high-lying states are suppressed by 
relaxation, and hence we have to take only the low-lying states into account if we consider 
an isolated quantum system as a model for the physical system. 
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The errors to the pumped current arise from the leakage current due to finite residual 
Josephson energy Ej'^'^ and from the non-adiabaticity of the temporal evolution which we 
consider here. Since the pumped current is determined by 2en™^^/T, the optimal way to 
pump is to choose the integer n^^^ to be as large as possible to decrease the contribution to 
the period T from opening and closing the SQUIDs, see Fig.|21 Thus the maximum pumped 
current is determined by the maximum value for rig limited by the criteria of Eqs. (|48p 
and (|52p when one of the SQUIDs is fully open and the other one is closed. Since the energy 
difference Eg — Sg depends strongly on rig in the two state approximation valid for 5 <C 1, 
the estimate in Eq. PHj) does not hold. In contrast, we adjust rig for all t such that 

|(e;q|(aii^si)|g;q)| ^^ .^^. 



£e{t) - eg{t) 

where dtHsi = 2_E'c%(% — ""■)• Substitution of Eqs. (jHHjl and (jH^jl into Eq. (|HHjl results in a 
differential equation for the gate charge 

ng=^-^[{ng-l/2r + S% (54) 

which implies together with the initial condition rrg(O) = that 

rrg(t) = - + Stan [iK^t - arctan((5~72)] • (55) 

The above equation yields the period T = {n/2 — 26)/ Hc + 0{P) and the pumped current 

I^,^ = -2e^Jl + ^]+0{6'), (56) 

where Kg is to be chosen such that Hc <^ re/4, see Eq. (J52|) . Although a high damping rate 
is desirable in hght of the adiabaticity, it requires a strong coupling of the system to its 
environment. The strong coupling has disadvantages, since it may introduce noise to the 
control parameters of the Hamiltonian, and hence result in pumping errors. 

Equation (jSUj) shows that the pumped current increases with 6 = E^^^/Ec- This ob- 
servation motivates us to work in the limit 5^1 for which the phase states are accurate 
eigenstates of the Hamiltonian in Eq. (Q). Since the ground state of the system is localized at 
the minimum of the potential, we may approximate the Hamiltonian Hs\ with a Hamiltonian 
of a perturbed harmonic oscillator as 

4ho = -Ecdl + ^{<P'f - Ej2 + Hp, (57) 

where 0' = — if/2, the perturbation Hamiltonian assumes the form Hp = —Ej2[cos{(f)') — 
1 + (0')^/2], and we have assumed for simplicity that Eji = 0. Furthermore, we have 
neglected the contribution from rrg since it can be taken into account in the unperturbed 
case by multiplying the wave function by e*"s<!^ which does not introduce energy shifts or 
changes into the matrix elements {k;q\ {dn //^gi) | r/i; q) . By solving the eigenstates of the 
unperturbed Hamiltonian Hpho — Hp and applying perturbation theory up to any order 
beyond the first it can be shown that 

^^V^ = v^"^ - 2-'(2m2 + 2m) + 0{6-'/'), (58) 

Ec 
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and 

(m;q|(g„^gsi)|0;q) ^ f i2^l\b^l^ - 2-^^-^!^) + 0{5-^/^) , for m = 1 
Ec \ 0{5-^/^) , for m ^ 1 



(59) 



In fact the eigenstates of the unperturbed Hamihonian are also the eigenstates of the per- 
turbed Hamiltonian up to (9(5"^/^), which justifies our assumption that n^ may be neglected 
also in the perturbed case. We showed in the above discussion that the assumptions leading 
to the bound in Eq. (jiH|) are valid here. Substitution of Eqs. (j^ and (j^ into Eq. (P5|) for 
171 = 1 yields 

\C,{t)\ < ^ [21/^5-3/^ + 3 X 2-13/45-5/4 ^ 0(^-7/4)1 _ ^ (gQ) 

Eq 

Hence, the pumped current in this limit Ej^'^^ ^ Eq is obtained from 

/p,„ = eh, = ^ [2-'/'6'/^ - 3 X 2-''/^6'/' + 0{6-'/')] , (61) 

ft 

where we choose A <^ 1 to assure that the system is accurately in its ground state. The 
leading contribution to the current may be expressed as 

2-'/'e\E%'El/' 
/pho ~ ^^^' (^2) 

where the Josephson energy Ej2 = E^^^- To justify the accuracy of above approximations, 
we show in Fig. El the estimates for the pumped current obtained from Eqs. (jU^ and (pT|) 
compared with a numerical solution, for which the spectrum of the Hamiltonian H^i was 
calculated, and the gate charge n, was ramped up employing Eq. (jIH)) and the condition 
|C*m(^)| < A for all m. On the other hand, the inset in Fig. El shows that dependence of the 
energy difference Si — £0 on Ug, and hence time, weakens with increasing 6, which renders 
Eq. (HHj) to be valid. Thus we conclude that Eq. (j^ yields an accurate analytic expression 
for the pumped current in a Cooper pair sluice for Ef^^^ ^ Eq- 

The calculation of the detailed dependence of the pumping errors 6Ip on the parameter A 
is beyond the scope of this paper. Instead, we will employ an approximation 

^ ^ A^. (63) 

-'p 

Thus we only give order of magnitude estimates for the maximum pumped current with a 
given error rate. 

Figure El shows that to maximize the pumped current, we should ultimately maxi- 
mize E^^^- Apart from choosing the material of the superconductor such that its energy 
gap is as large as possible and using very transparent tunnel junctions, the only way to 
increase the maximum Josephson energy of the SQUIDs is to increase the surface area of 
the tunnel junctions. However, the junction capacitance increases with the surface area, and 
hence the Coulomb energy Eq decreases. On the other hand, the Coulomb energy should 
be large compared with the temperature for the finite temperature effects to be minimal. 
Thus the attainable base temperature limits E^^^- For standard experimental parameters, 
i.e., relatively transparent aluminum junctions, Ec/ks = 1 K, and for A = 0.01 we obtain 
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FIG. 3: Pumped current as a function of the maximum Josephson energy of the SQUIDs E 
obtained from Eq. (|E^ (dotted hue), Eq. (|61() (dashed), and full numerical solution of the system 
(solid line). The inset shows the maximum relative variation of the energy difference e\ — eo with 
respect to rig as a function of the maximum Josephson energy of the SQUIDs. 



currents of the order of a nanoampere. The pumping errors are not critical for observation 
of the Berry phase, and hence the parameter A can be chosen to be, e.g., 0.1 which implies 
lower than one percent population in the first excited state and a pumped current of the 
order of 10 nA. On the other hand, if the sluice is used to define the current standard, 
relative pumping errors should be less than 10~^ implying A < 3.1 x 10~^. In this case, 
the standard parameters still yield Jpho to be of the order of 10 pA. For superconducting 
niobium, the gap in the quasiparticle energy spectrum is an order of magnitude higher than 
for aluminum leading to about six times larger pumped currents than for aluminum devices. 

V. MEASUREMENT SCHEME FOR THE PUMPED CURRENT 

To observe the pumped current of a phase biased Cooper pair sluice, we install a third 
tunnel junction and a DC current source with current Jdc in parallel with the sluice, as shown 
in Fig. I3J The applied current Jdc is assumed to be smaller than the critical current of the 
system, and hence there is no voltage over the third junction. In addition, the capacitance of 
the third junction is assumed to be large enough such that its charging energy is negligible. 
Thus the third junction is essentially in such an eigenstate of its phase operator ip-^, that the 
free energy of the system is minimized. We write the free energy of the circuit shown in 
Fig. H as 

G = Hsl- E]3 cos V33 — Vs, (64) 

where the second and the third term correspond to the Josephson energy of the third junction 
and to the work done by the current source, respectively. Since the system is in its ground 
state in the adiabatic temporal evolution, the free energy must achieve its minimum at (^3, 
leading to 

d^^g = Ejs sin ips + —(/si - Jdc) = 0, (65) 

where Jgi = {Ejisimpi + Ej2siia{p2)/2 denotes the total average supercurrent through the 
sluice, and we have used the relation 27r$/$o = V^i + ¥^2 — V^s- In Eq- dHHl), we did not, 
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however, include the effective current Jp = 2erif^^ /T through the third junction arising 
from the pumping. This current has to be subtracted from the apphed DC current and thus 
we obtain a relation 

/ := Jdc - -^p = — z — sin Lp^ + hi- (66) 

Equation ()66p shows that the pumping of Cooper pairs alters the current through the system, 
which can be measured as a shift in the critical current of the device. 




FIG. 4: Circuit diagram for the measurement of the pumped current through the Cooper pair 
sluice. The third tunnel junction and a DC current source with current /dc is installed in parallel 
with the sluice shown in Fig. ^ The phase difference of the third junction is denoted by ^p^. 



In the case of ideal SQUIDs, for which e = Ef^ /Ef^^ = 0, it follows from Hamiltonian 
of Eq. (j37|) that the pumped current is independent of the total phase difference ip. Thus 
we obtain from Eq. (j66j) 

/= smypg, (67) 

which is equivalent to the case where the sluice is replaced by an ideal current source. Since 
the pumping shifts the effective current I through the third junction by /p, the pumping can 
be observed by switching measurements, i.e., by applying DC current pulses of duration r 
and different magnitudes /dc while monitoring the voltage across the device. If the device 
switches to the normal conducting state, a non-zero voltage is observed. 
We estimate the switching probability in the zero temperature limit as 



Ps, 



-TrMQT 



(68) 



The transition rate for the macroscopic quantum tunneling in a weakly dissipative environ- 
ment is given by 

/,!. I 36/V7 

(69) 



LOr 



Fmqt = l2VQi-^JMJ/{nuj^)e 



27r 



where Af/ ^ \E^:>,[2{1 - I / Icz)f/^ and 



UOr 



8Ej3/;c3[2(l - I/Ic3)V/yh m denote the 



barrier height and the plasma frequency of the 3rd junction, respectively. Figure El shows 



16 



the switching probabihty given by Eq. ()68|) as a function of the DC current Jdc with forward 
and backward pumping for pumped currents of 10 nA and 1 nA which were found to be 
experimentally feasible in Sec. IIVI In Fig. Efa), we have assumed that /cs = 0.1 /xA and 
Cj3 = 200 fF. Since the sluice is also assumed to be optimized for high Josephson energy, the 
additional assumption of the parameters of the third junction implies that for example in 
standard shadow evaporation combined with electron beam lithography, three-angle evapo- 
ration is to be employed to obtain a different transparency of the third junction as compared 
with the junctions used in the sluice. An alternative to the three-angle evaporation is to 
install an additional capacitor in parallel with the third junction, as realized in Refs. (l5l.ll6l|. 
Although these techniques in the fabrication of the sample are very feasible, we have plotted 
in Fig.El^b) the switching probability without them, i.e., for /c3 = 2 /iA and Cj3 = 15 fF. To 
detect the pumping for the parameters of Fig. Efb), the relative visibility in the variation of 
the applied DC current should be of the order of 0.1% for /p = 1 nA and 1% for Jp = 10 nA. 
We note that in our current experimental set up for other switching measurements than the 
proposed one, we have detected smaller than 0.1% shifts in the current histograms similar 
to the ones shown in Fig. |3| Thus an accurate measurement of the pumped current for the 
parameters of Fig. El^a) should be achievable. 



(a) 



(b) 



1 




/ 




fff 


0.8 






: 0.8 


1 1 


0.6 






0.6 

s 

a,- 


II 


0.4 






0.4 


: ; 


0.2 


1 ; 




0.2 


w 





; JJ 




Q 


^■'' 


60 80 


1400 1450 1500 1550 




^DC[I1 


A] 




/dc [nA] 



FIG. 5: The probability of the system of Fig. ^ to switch to the normal conducting state as a 
function of current Jdc with 10 nA backward pumping (dashed Une), 1 nA backward pumping 
(dash-dotted line), 1 nA forward pumping (soHd hne), and 10 nA forward pumping (dotted line). 
The parameters of the third junction are chosen to be (a) /cs = 0.1 \ik. and Cjs = 200 fF, and (b) 
/c3 = 2 /iA and Cjs = 15 fF. The duration of the DC current pulse is r = 100 /is. 



The correction to the shift of the critical current due to finite residual Josephson en- 
ergy E^^ is of the order of max(/ci, /c2)e which is small compared with the pumped current /p 
for typical experimental parameters. Thus we conclude that the proposed measurement 
scheme provides a convenient experimental arrangement for observation of Berry phase in 
superconducting circuits. 

Furthermore, we note that Fig. 0] strongly resembles the circuit used in the Saclay experi- 
ments on the artificial two-state atom quantronium [lall^- The sluice in Fig. |3] corresponds 
to the quantronium and the other part of the circuit is introduced for measuring the su- 
percurrent through it. Since the supercurrent depends on the state of the quantronium, 
the switching of the superconducting system into normal state can be employed to observe 
the state of the quantronium. Thus the supercurrent cannot be neglected as we did in the 
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derivation of Eq. ()67|). and the exact calculation of the switching becomes more difficult 



than in our case. For further details, see Refs. |22, |2 



VI. CONCLUSIONS 



In this paper, we analyzed a phase biased Cooper pair sluice, for which the pumped 
charge is closely related to the Berry phase accumulated in one pumping cycle. We showed 
that the adiabaticity requirement of the sluice does not restrict the pumped current to 
be undetectably small if the sluice is optimized according to the results derived. In fact, 
pumped currents of tens of nanoamperes might be achievable. Moreover, a lower bound for 
the maximum pumped current allowed by the adiabaticity criterion for relative pumping 
error of 10~^ and standard experimental parameters turned out to be of the order of 10 pA. 
In Ref. jl3|, a pumped current of 100 pA was obtained with the error rate 10~^ using a 
dynamical approach, which suggests that the lower bound for the maximum pumped current 
obtained in this paper can be increased by a more detailed derivation. Thus the sluice is a 
potential candidate for a metrological current standard. 

We presented a detailed scenario for observing the pumped current, and showed that 
the visibility of the current in the proposed measurement is expected to be high. In future 
work, we aim to experimentally observe the Berry phase employing the scheme given in this 
paper. Although we have neglected sources of error, e.g., decoherence in our analysis, we 
expect that they are not dominant in the experimental realization 26|. In fact, our goal 



is only to combine two already realized techniques: the Cooper pair pumping [l^ and the 
switching measurement with a large Josephson junction as an ammeter. Thus the scheme 
for observing the Berry phase in superconducting circuits presented here is considered to be 
very promising. 
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